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Abstract 

We investigate the viability of observing charged Higgs bosons (H^) produced 
in association with W bosons at the CERN Large Hadron Collider, using the lep- 
tonic decay H + — * t + v t and hadronic W decay, within different scenarios of the 
Minimal Supersymmetric Standard Model (MSSM) with both real and complex 
parameters. Performing a parton level study we show how the irreducible Stan- 
dard Model background from W + 2 jets can be controlled by applying appropriate 
cuts and find that the size of a possible signal depends on the cuts needed to sup- 
press QCD backgrounds and misidentifications. In the standard maximal mixing 
scenario of the MSSM we find a viable signal for large tan j3 and intermediate 
masses (~ mj) when using softer cuts ($±,p± T . iet > 50 GeV) whereas for harder cuts 
(tf±,P±T- Jct > 100 GeV) we only find a viable signal for very large tan/3 (> 50). We 
have also investigated a special class of MSSM scenarios with large mass-splittings 
among the heavy Higgs bosons where the cross-section can be resonantly enhanced 
by factors up to one hundred, with a strong dependence on the CP-violating phases. 
Even so we find that the signal after cuts remains small except for small masses 
(< m t ) when using the softer cuts. Finally, in all the scenarios we have investigated 
we have only found small CP-asymmetries. 



I Introduction 

The quest for understanding electroweak symmetry breaking and generation of elementary 
particle masses is one of the driving forces behind the upcoming experiments at the CERN 
Large Hadron Collider (LHC). In the Standard Model (SM) the electroweak symmetry 
is broken by introducing one Higgs doublet whereas in the Minimal Supersymmetric 
Standard Model (MSSM), which is the main focus of this paper, two Higgs doublets 
are needed. More specifically, the MSSM is a two Higgs Doublet Model (2HDM) of type 

II where one of the doublets gives mass to all the up-type fermions and the other to all the 
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down-type fermions. After electroweak symmetry breaking the Higgs sector in the MSSM 
consists of three neutral Higgs bosons and one charged. The charged Higgs boson (H^) 
is of special interest since there are no charged scalars in the SM and thus its discovery 
would constitute an indisputable proof of physics beyond the Standard Model. Therefore, 
the hunt for charged Higgs bosons will play a central role in the search for new physics at 
the LHC experiments. 

Currently the best model-independent direct limit on the charged Higgs boson mass is 
from the LEP experiments, m H ± > 78.6 GeV (95% CL) [1] (assuming only the decays 
H + — > cs and H + — > t + v t ). In addition there are stronger direct limits from the Fer- 
milab Tevatron, but only for very large or very small tan/3 (tan/3 > 100 and tan/3 < 1, 
respectively) where tan (3 denotes the ratio of the vacuum expectation values of the two 
Higgs doublets, as well as indirect limits mainly from B-decays although the latter are 
quite model-dependent. We refer to the PDG [2] for details. 

The main production mode of charged Higgs bosons at hadron colliders is in association 
with top quarks through the gb — > H~t and gg — > H~tb processes [3-6] with the former 
one being dominant for heavy charged Higgs bosons m H ± > m t and the latter one for light 
ones m H ± < m t —m b . The possibilities of detecting charged Higgs bosons in these channels 
have been studied extensively (see for example [7-13] and also more specifically for the 
two LHC experiments [14-16]). In summary, the leptonic decays H + — > t + v t seem to be 
most promising for detecting light charged Higgs bosons at any value of tan (3 and heavy 
charged Higgs bosons for large tan (3 whereas the hadronic decays H + — > tb may be useful 
both for large and small tan (3 above threshold. However, the transition region m H ± ~ m t 
has so far been difficult to cover and the same is true for the intermediate tan/3 ~ 
^rrit/mb region, except in special supersymmetric scenarios with decays into charginos 
and neutralinos [17-19]. The possibility of using the gg — > H~tb process to improve the 
discovery potential in the transition region has been studied in [20,21]. Recently a new 
method for matching the differential cross-sections for the two production modes has been 
developed [22] resulting in a significantly improved discovery potential in the transition 
region [23]. 

Given the problems with exploiting the production of charged Higgs bosons in association 
with top-quarks it is also important to investigate other production modes. One primary 
example is the production in association with W bosons [24-33]. In addition to being a 
complement to the main production mode, especially in the transition region, this channel 
may also give additional information on the Higgs sector. So far, H ± W Zf production has 
mainly attracted theoretical interest with limited attention to the experimental viability, 
the only exception being [26] which came to the conclusion that in the MSSM the signature 
W ± H T — > W ± tb — > W ± W Zf bb cannot be used due to the irreducible background from 
tt production. On the contrary, in this paper we will show that the situation may be 
drastically improved, at least for large tan (3, by instead looking at leptonic decays H + — > 
t + v t together with W — > 2 jets. Alternatively, going to a more general 2HDM one can 
get a very large enhancement of the production cross-section compared to the MSSM 
through resonant enhancement [33]. The associated production of charged Higgs bosons 
with CP-odd scalars A and subsequent leptonic H + — > t + v t decays have been analyzed 
in [34]. 

Different types of higher order corrections to H ± W T production that have been studied 
typically give effects of the order 10-20%. This includes the supersymmetric electroweak 



corrections [30], the next-to- leading order (NLO) QCD corrections [31], and the supersym- 
metric QCD corrections [32]. Of special importance here are the NLO QCD corrections 
since they show that one has to use running masses for the Yukawa couplings, which in 
turn has a large impact on the magnitude of the cross-section as will be discussed later. In 
addition to process specific higher order corrections there are also higher order effects from 
the bottom/sbottom sector in supersymmetric theories which are large for large tan /3 and 
therefore have to be resummed to all orders in perturbation theory [35,36]. This results 
in additional terms proportional to Amt, in an effective Lagrangian [36] describing the 
couplings of charged Higgs bosons to top and bottom quarks. Finally we expect that a 
variation of the SUSY scenario, especially of /x, will lead to a variation of the cross-section 
of 0(10%) [37]. Since we will only be studying the H ± W T production at parton level, 
although with appropriate smearing of momenta, we have chosen not to take into account 
any of these corrections in the production, except for the use of a running 6-quark mass 
in the Yukawa coupling. 

The study of pp — > H ± W T also offers the possibility to explore effects of CP violation. 
In the MSSM with complex parameters one gets new sources of CP violation beyond 
the CKM matrix which can give rise to explicit CP violation also in the Higgs sector 
through loops of supersymmetric particles [38,39]. For example, such phases may give 
rise to differences in the — > W + H~ and — > W~H + cross-sections which 
have been analyzed in a more general CP-violating 2HDM [40] although not taking into 
account non-diagonal Higgs propagator effects [41,42]. Another example is the CP-odd 
rate asymmetry induced by loops of SUSY particles in the decay H + — > tb [43] that has 
been analyzed for the LHC in [44]. 

The outline of this paper is as follows: We start out by giving the differential cross-section 
for the dominant production mode of charged Higgs bosons in association with W bosons 
for large tan/5 and study the dependence on m H ± and tan/? at the LHC. In section 3, 
we then specialize to the decays H + — > t + u t , W~ — > 2 jets and compare signal and 
background giving special attention to the difficult transition region where m#± ~ m t 
and the leptonic decay mode of the charged Higgs boson is dominant. We show how the 
irreducible SM background for the resulting signature rp 1 ^ + 2 jets can be substantially 
reduced by appropriate cuts. In section 4, we then give the results of our study in the 
MSSM both with real and complex parameters. In the latter case we investigate the 
possible dependence of the signal on some of the CP-violating phases of the MSSM. 
We also study a class of special scenarios within the MSSM where the cross-section is 
resonantly enhanced. Finally, section 5 contains our conclusions. 

2 H±W^ production at the LHC 

In this section we give the differential cross-section for the dominant production mode 
of charged Higgs bosons in association with W bosons at large tan/5 in a general type 
II 2HDM and then specialize to the case of the maximal mixing scenario in the MSSM 
giving the dependence of the total cross-section on m H ± and tan/3 at the LHC. 




2.1 Differential cross-section 



In a general 2HDM of type II, the dominant production mechanisms for a charged Higgs 
boson in association with a W boson at hadron colliders are bb annihilation 

bb -> H^W* (1) 

at tree-level and gluon fusion 

99 - (2) 

at one-loop- level [24,25]. In this study we focus on the parameter region with intermediate 
H ± masses (~ m t ) and large tan/5, where the decay H ± — > tv t has a large branching 
ratio and where the bb annihilation dominates, and thus we do not need to consider the 
contribution from gluon fusion. Note, that this is true also in MSSM scenarios with large 
tan/3, light squarks and large mixing in the squark sector because the possible strong 
enhancement of the gluon fusion cross-section through squark loops occurs only for small 
tan/5 < 6-7 [29]. 

The Feynman diagrams for bb — > H ± W Zf are shown in figure 1. The bb annihilation 
proceeds either via s-channel exchange of one of the three neutral Higgs bosons in the 
2HDM or t-channel exchange of the top quark. 

In the real 2HDM there are two CP-even, h° and H°, and one CP-odd, A , neutral Higgs 
bosons, whereas in the general complex case all three neutral Higgs bosons mix and form 
three mass eigenstates Hi, i = 1, . . .3. The Lagrangians relevant for the bb — > H + W~ 
process are, following the conventions of [40,45,46], 

C Hi H±w* = 9 -Y.9H i H-w + {H i id^H-)W + ^ + h.c, (3) 



CHlib = ~ 2m w cos 13 J2 Hib ^ 9H ^ ~ 1 lm (9H-bbh5]b, (4) 

i 

C H ±tb = -J— H + i [m t cot (3 hi --f 5 )+m b tanp hi + l5 )}b + h.c. (5) 
V 2mw 1 1 



with the couplings 



9HiH-W+ 

9Uilb 



= 9* Hi h+w- = 2i cos(3- 
= On + i 3i sin (3 , 



On sin (3 + i 3i , 



(6) 



where Oji is the Higgs mixing matrix. In the real 2HDM, Hi = {h°, H°, A } and the 
mixing matrix has the simple form 



which gives purely real couplings for h°, H° and imaginary ones for A . 

In the general CP-violating 2HDM all entries in the mixing matrix can be non-zero. The 
resulting mixing between CP-even and CP-odd Higgs states leads to complex couplings in 
general, which together with the propagators gives differences between the cross-sections 
for bb — > H + W~ and bb — > H~W + in analogy to the — > W ± H T processes analyzed 
in [40]. However, it has to be kept in mind that the formalism we are using is simplified 
and only valid if at most one of the neutral Higgs propagators can be resonant. Otherwise 
one has to use the complete description [41,42], which also takes into account non-diagonal 
propagator effects that arise from the mixing of two different Higgs bosons through higher 
order loops. Using the simplified formalism, the differential cross-sections for the two 
processes are [25,40]: 
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where s and t are the ordinary Mandelstam variables of the hard process and 

X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zx) , (10) 
2 _ A(g,m 2 y ,m 2 ^ ± )sin 2 ^ 

with 9 being the polar angle in the 2-^2 cms and 



In the remainder of this paper we focus on H^WT production in the MSSM, with real 
and complex parameters, which is a special case of a type II 2HDM. The Higgs masses, 
widths, branching ratios and the mixing matrix of the neutral Higgs bosons in the MSSM 
can be calculated with programs such as FEYNHlGGS [47-50] and CPSUPERH [46], which 
include higher order corrections up to leading two-loop contributions, especially the Amj 
corrections. 

As already alluded to, H ± W Zf production could potentially be used to determine the 
amount of CP violation in the MSSM. In order to isolate the possible effects of CP 
violation, a CP-odd rate asymmetry 



iCP 



a(pp -> H+W~) - a(pp -> H~W + ) 



a(pp -> H+W-) + a(pp -> H-W+) 
can be defined, where only interference terms contribute to the numerator 

a(pp -> H + W~) - a(pp -> H~W + ) oc ^-{bb -> H + W~) - ^-{bb -> H~W + ) = 
= 12tts | c^/? 2^^9H i H-w+9H j H-w+m{S Hi S H .)BB{g H . ib g Hj 



(13) 



+ 



— b — [m 2 w m 2 H± - spl - t 2 ] 22lm(g HiH - w+ g Hilb )lm(S Hl ) )■ (14) 

m t ) cos p . j 



From this expression it is clear that the conditions for obtaining large CP asymmetries 
are quite subtle. Starting with the first term, which comes from the interference of two 
different s-channel amplitudes, we see that in this case large asymmetries can be obtained 
only if the phases in the propagators make Im(S'^S'l^) large and at the same time the 
Higgs mixing matrix Oij leads to large factors from the couplings Yvci(gH i H-w+9*H :j H-w+) 
and Ke(g H .^ b g^j r b ). Typically this means that the maximal asymmetry with one resonant 
propagator is of order T H ./m H ., modulo the coupling factors. In the more general case 
with several resonant propagators the formalism we are using is not sufficient as already 
noted. Looking at the second term in Eq. (14), which originates in the interference of 
one s-channel amplitude with the t-channel one, this can give an asymmetry of the order 
THi/iTLHi if the s-channel amplitude is resonant, again modulo the coupling factor. How- 
ever, in this case the coupling factor will make the asymmetry suppressed with 1 / tan (3 
for large tan (3 compared to the asymmetry arising from the interference of two s-channel 
amplitudes. 



2.2 Cross-section calculation 

We have implemented [51] the two processes bb — > H + W~ and bb — > H~W + as separate 
external processes to Pythia [52, 53] in order to be able to analyze the rate asymmetry 
Aq-p. In principle, the implementation in Pythia makes a generation of the complete 
final state possible, but for this first study we have chosen to stay on leading order parton 
level. 

We use running b and t masses in the Yukawa couplings and a running electroweak coupling 
aEw, all evaluated at the renormalisation scale /j,r = rriH± + mw using the leading order 



MSSM Parameters. All masses in GeV. 



Scenario 


A t 


A b 


M 2 


m- g 


m h 


m h 




m b 2 


Maximal mixing (m™ ax ) 


2000 


2000 


200 


800 


820 


1177 


996 


1012 


Less mixing 


1000 


1000 


200 


800 


922 


1099 


993 


1012 


No mixing 








200 


800 


1014 


1015 


991 


1014 



Table 1: MSSM parameters for the maximal mixing scenario, m™ ax , the less mixing 
scenario and the no mixing scenario in addition to tan (5 — 50, /i — 200 GeV , and 
Msusy = 1000 GeV, as well as the resulting 3rd generation squark masses. 

formulas implemented in Pythia. Numerically, for = 175(400) GeV the masses are 
rrib = 2.7(2.6) and m t = 160(153) GeV. This leads to a reduction of the cross-section 
compared to the naive tree-level without running to about a third, which is in better 
agreement with the NLO calculation [31]. As factorization scale we use the average mass 
of the H ± and the W, Hf = ^( m H ± + n^w), again for better agreement with NLO 
calculations. All other parameters are left at their default values in Pythia, for example, 
we use the CTEQ 5L [54] parton densities. We note that the main uncertainty in the 
6-quark parton density is not so much a question of which set that is used, but rather 
the precise value of the factorisation scale. Varying the scale with a factor 2 up or down 
the cross-section changes with ±20%. This should be compared with the about ±5% 
uncertainty in the b density given for example by CTEQ65E [55] for relevant values of x 
and Q 2 . 1 In addition the widths of the H* 1 and W bosons are also included in the same 
way as in standard PYTHIA. In other words the m#± and mw masses vary according 
to Breit-Wigner distributions with varying widths meaning that for each mass the decay 
widths are recalculated based on the actually open decay channels. 

For the calculation of the MSSM scenario and the corresponding Higgs masses, the Higgs 
mixing matrix Oji and the branching ratios of H ± we use FeynHiggs 2.2.10 [49]. From 
Oji we then calculate the couplings gUiH-w+ an d g^U according to Eq. (6). 

Figure 2 shows the resulting total cross-section a(pp — > H ± W T ) as well as the H ± branch- 
ing ratio into tv t as functions of the H ± mass and tan/? at the LHC. Here, and in the 
following unless otherwise noted, we have used a maximal mixing scenario, defined in ta- 
ble l 2 . At the same time, we have found that the cross-section depends very little on the 
mixing in the third generation squark sector. For example, for a mass of m#± = 175 GeV 
and tan (3 = 50 the cross-section decreases from 260 fb to 257 fb in the no mixing scenario 
of table 1. 

Apart from the decreased cross-section due to the use of a running mass instead of the 
pole mass, our results also differ from, for example, the ones by [25] at large Higgs masses 
due to off-shell resonant enhancement. When the width of the H ± is large, the mass of 

1 For this estimate we have used the online tool for pdf plotting and calculation at http://durpdg. 
dur . ac . uk/hepdata/pdf 3 . html. 

2 Here tan/3 = Vijv\ denotes the ratio of the vacuum expectation values of the Higgs fields and /j, 
is the Higgs/higgsino mass parameter. Msusy defines a common value of all squark and slepton soft 
SUSY breaking mass parameters Mq, Mfj, M 3 D , M£, MJL, where g = 1,2,3 is the generation index. 
X t — A t — \i* cot (3 and Xt, = At, — fi* tan (3 describe the mixing in the third generation squark sector with 
the trilinear scalar couplings A t and A^. In the gaugino sector M 2 denotes the SU(2) soft SUSY breaking 
mass parameter, nig is the mass of the gluinos and the U(l) soft SUSY breaking mass parameter Mi is 
fixed by the GUT relation M 1 /M 2 = 5/3 tan 2 6 W . 
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Figure 2: The left plots show the (bb contribution to the) cross-section a(pp — * H ± W Zf ) 
at the LHC and the branching ratio 5i?(/J ± — ► rv T ) as a function of m H ± for tan /3 — 50 
(solid line) and tan/5 = 10 (dashed line). The right plots show the cross-section and 
branching ratio as a function of tan (3 for m H ± = 175 GeV (solid line) and m H ± = 400 GeV 
(dashed line). 
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Figure 3: Distribution of the invariant mass of the decay products (in our case rv T ) 
from the H ± showing the off-shell enhancement of the cross-section for m H ± = 600 GeV, 
resulting in m,H 123 = 136,594,594 GeV. In this case the large width Y H ±(m H ±) = 20.3 
GeV means that the mass can fluctuate down to m TVT + m w < m H . which results in the 
broad peak seen to the left. 



the decay products (in our case ru T ) can fluctuate down such that m TVr + mw < rn^ 
giving rise to s-channel resonant production. This can be seen from the big broad peak 
in the m TVr invariant mass spectrum shown in figure 3, which has been obtained for the 
following masses: m H ± = 600 GeV, = 136 GeV, m# 2 = 594 GeV, and m# 3 = 594 
GeV. In turn this leads to an enhancement of the cross-section for large Higgs masses. 
The same effect can also be seen in the study by [26]. However, it is not clear to what 
extent this increase in the cross-section leads to a larger signal in the end since the much 
wider peak could potentially be harder to see. This is especially true for the leptonic 
decays we are considering since here we can at best reconstruct the sum of the transverse 
momenta carried by neutrinos. 



3 Signal selection 

In this section we explain the final signature analyzed in this paper, discuss background 
processes and define the cuts necessary to enhance the signal-to-background ratio. 

The signal selection is based on the intermediate mass =175 GeV for tan/5 = 50 in 
the maximal mixing scenario given in table 1. Results based on this selection for other 
values of the mass and tan (5 as well as in other scenarios will be given in the next section. 



3.1 Signature 

As already mentioned, in this study we focus on associated H ± W T production with 
subsequent leptonic decays of the charged Higgs, H + — > t + u t and hadronic decays of the 
W boson, W~ — > qq', q — u,c, q' — d, s. The decays W + — > t + v t and H — > cs, cb 
have not been included since the branching ratios for the latter are negligible compared 



to H + — > t + v t for large tan/5. For simplicity we only consider hadronic decays of the 
r-lepton, r — > z/ r + hadrons. The resulting signature thus consists of two light jets (J), 
one hadronic r jet (rj e t) and missing transverse momentum (jzfjj carried away by the two 
neutrinos (one from the charged Higgs boson decay and one from the r-lepton decay): 

2j + Tjet + i>L ■ 

Due to the two neutrinos escaping detection it is of course not possible to reconstruct the 
invariant mass of the charged Higgs boson in this decay mode. However, from tf± and the 
transverse momentum of the r jet, p_i_ Tjct , the transverse mass 3 

m± = i y2p ±rjot ^[l-cos(A0)] (15) 

can be calculated, where Acf) is the azimuthal angle between p± Tjct and tf±. If there is a 
detectable charged Higgs boson it will show up as a peak in this distribution with the 
upper edge of the peak given by the mass of the charged Higgs boson. 

3.2 Background 

The dominant irreducible SM background for our signature 2j + r- jCt + arises from W + 
2 jets production which we have simulated with help of the package ALPGEN [56]. The 
program calculates the exact matrix elements on tree-level for the 2j + r + u T final state. 
In this way it includes not only W + 2 jets production as well as W pair production, but 
also contributions where the r lepton and the neutrino do not arise from the decay of a 
(virtual) W boson. The obtained cross- sect ions and distributions of momenta, invariant 
masses etc. have been cross-checked with Madgraph [57-59]. 

As a further precaution we have also simulated another irreducible background, namely 
WZ + 2 jets production with Z — > vv using ALPGEN, which typically has a larger missing 
transverse momentum. However, we have found that this process contributes less than 
3% to the overall background after appropriate cuts (table 2 below) have been applied. 
So even though the tails of the p± and m± distributions are slightly harder compared to 
the W + 2 jets background, it is still safe to neglect the WZ + 2 jets background for this 
study. 

The importance of using the complete matrix element for the signature we are considering, 
and not just the Wjj approximation, can be seen in figure 4. The figure shows the 
invariant mass m TUr of the r and v T system generated with ALPGEN compared to the W 
mass from Pythia in H^ + jet production, which has a normal Breit-Wigner distribution 
with varying width. Of special interest here is the tail for large invariant masses m TVr and 
we note that there is a factor of 2 difference between the two approaches in this tail. To 
illustrate the invariant mass range of relevance in our case we also show m TUT generated 
with ALPGEN after all event selection cuts (see table 2) have been applied. This also 
shows that the SM background can be considerably reduced by appropriate cuts, which 
will be discussed in the next subsection. 

3 Strictly speaking this is not the transverse mass since there are two neutrinos in the decay chain of 
the charged Higgs boson we are considering. Even so the characteristics of this mass is very similar to 
that of the true transverse mass. 




Figure 4: Invariant mass m TVr of the r and v T system generated with ALPGEN (solid line) 
compared to the W mass from Pythia in W + jet production (dotted line), which has a 
normal Breit-Wigner distribution with varying width. The distributions have been nor- 
malised to the same height. Also shown is the m TVr distribution generated with ALPGEN 
after applying the cuts defined in table 2 (dashed line). 



3.3 Cuts 



Our study is performed at parton level, without any parton showering or hadronisation. 
Instead the momenta of the jets are smeared as a first approximation to take these, as 
well as detector effects, into account. The only exception is that for both signal and 
background Tauola [60, 61] has been used to perform the decay of the r lepton into 
a hadronic jet plus a neutrino. The use of Tauola also makes it possible to take into 
account the differences in the decay characteristics of the r-lepton depending on whether it 
comes from a charged Higgs bosons with spin or a W-boson with spin 1. This difference 
can, at least in principle, also be used to discriminate against backgrounds [9,11,62,63] 
as discussed below. 

We distinguish between the two jets by calling them hard (with momentum phj) and 
soft ip S j) according to the larger and smaller value of their transverse momentum p±, 
respectively. We have done a Gaussian smearing of the measurable momenta p T . ct , ptj and 
p S j and then calculated^. The smearing is done preserving the direction of the momenta 

and the width of the smearing is defined as a [GeV] = y ^gev] + & with a = 0.60 and 
6 = 0.02. 

Our set of basic cuts, see table 2, is defined according to the coverage and resolution of a 
realistic detector. First we require that the pseudorapidity of the r- ]Ct and the light jets are 
in the range —2.5 < rj < 2.5. A further cut requires a minimal distance ARjj > 0.4 and 
AR Tjct j > 0.5 between two respective jets to allow their separation. Finally we require a 
minimum p± > 20 GeV for all jets. Figure 5 shows the resulting rrijj, m± distributions 
and the p± distributions for all jets and the missing transverse momentum for the signal 
and background after these basic cuts. 

Based on the shape of the signal and background in figure 5 we define further cuts. These 
cuts are used to reduce the background while keeping as much signal as possible. It can 
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Figure 5: Top row: rrijj and m± distributions for signal (solid), and background (dashed). 
Middle row: p± distributions for Tj Ct and jf. Bottom row: p± distributions for the hard 
and soft light jet. In all plots the basic cuts of table 2 have been applied. 



be seen that a cut 70 GeV < rrijj < 90 GeV will be very efficient since the signal is 
peaked around the W mass whereas the background is almost flat. A cut m± > 100 GeV 
will also remove a large part of the background. Increasing this cut leads to smaller 
S/y/B for m H ± = 175 GeV, yet might be useful for higher m H ±. The signal peak has 
a soft upper edge, so for Higgs masses below 125 GeV it would be very hard to see the 
signal since the background increases dramatically in this region. To further reduce the 
background the cuts p±hj > 50 GeV, p± S j > 25 GeV are defined. In addition, we apply 
the cuts p± T j Ct > 50 GeV and^j_ > 50 GeV, similarly to [34], in order to reduce the QCD 
background and the effects of detector misidentifications although we have not simulated 
these effects explicitly. Since these cuts may be too soft we will also show results for 
the harder cuts p±r iat > 100 GeV and tf± > 100 GeV, which have been used for example 
in [15]. 

As mentioned above the r-jet properties of the signal and background are also different 
due to the difference in spin of the and W bosons. The relevant measure is the 
ratio, R, between the transverse momentum of the leading charged ir, p± n , and the total 
transverse momentum of the r jet, 

R = ^ 

We assume that p± w is measured in the tracker independently of the transverse momentum 
of the r-jet and in order to take into account the tracker performance we apply Gaussian 
smearing on l/p± n with 

i r~ 22 2 

a(— MTeV- 1 ] = JO.522 + , , (17) 

P±tt V (rl^ [GeV]) 2 sin 6^ 

where 9 n is the polar angle of the ir [64]. 

Having applied all the cuts outlined above, the resulting distribution of R is shown in 
figure 6 both when keeping all hadronic r-decays as well as only selecting 1-prong decays. 
As can be seen from the figure, the differences between the signal and background is clear 
in both cases. Even so, we have found that applying cuts on R in either of the two cases 
does not lead to an improved overall significance S/ \/B compared to the case when all 
hadronic r-decays are kept and no cut on R is applied. (When selecting only 1-prong 
events a cut R > 0.8 turns out to be advantageous, but the overall significance still drops 
compared to keeping all hadronic r-decays.) At the same time it should be kept in mind 
that these conclusions may change when making a complete simulation of the final state 
including parton showers and hadronisation as well as a full detector simulation. It may 
also be possible to enhance the use of r-polarization by considering the 3-prong decays 
separately and looking at the sum of the transverse momenta of the like-sign pair instead 
of that of the leading charged pion [63]. 

In table 2 all basic and additional cuts are summarized and in table 3 it is shown how 
applying the different additional cuts one after the other affects the integrated cross- 
section. The signal is given for the m™ ax scenario in table 1 with a charged Higgs mass 
of 175 as well as 400 GeV and tan/3 = 50. We have used an integrated luminosity of 
300 fb _1 and a r detection efficiency of 30% to calculate S/y/B. Finally we note that, 
when the charged Higgs mass is large it can be advantageous to use the harder set of cuts 
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Figure 6: Comparison of R = pA_ 7T /PA.T jct between signal (for m™ ax with m H ± = 175 GeV 
and tan (3 = 50) and background. The left plot is for all events and the right for 1-prong 
events. The dashed (dotted) curve denotes signal (background). 



Basic cuts Additional cuts [all in GeV] 
/'..,,. .- 50. /f - 50 
70 < rrijj < 90 
m ± > 100 

P±hj > 50, p ±sj > 25 



Table 2: Set of basic cuts which define a signal region that corresponds to the sensitive 
region of a real detector and additional cuts which reduce the background according to 
the distributions in figure 5 as well as suppress QCD background and detector misidenti- 
fications. 
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560000 


55 


4900 


0.7 


3.3 


300 0.04 


P±r ict > 50 , ji ± > 50 


22000 


25 


2200 


1.6 


2.7 


240 0.2 


70 < rrijj < 90 


1700 


21 


1900 


5 


2.2 


200 0.5 


m± > 100 


77 


15 


1400 


16 


2.1 


190 2.3 


P±hj > 50 , p ±sj > 25 


28 


9.3 


840 


17 


1.5 


135 2.6 



Table 3: The effect of the different cuts on the integrated cross-section for background 
((7b) and signal (a s ) in the m™ ax scenario with m#± = 175 and 400 GeV for tan/? = 50 
(see table 1) as well as the number of signal events S and the significance S/ \f~B assuming 
an integrated luminosity of 300 fb -1 and a r detection efficiency of 30% . 

P±r jct > 100 GeV and|f_i_ > 100 GeV as illustrated below. In order to get a rough estimate 
of how important higher order corrections can be on the resulting significances we consider 
the following worst case scenario. Assuming the total uncertainty in the signal as well as 
background to be a factor 1.5, then the significances in table 3 would in the worst case 
be reduced with a factor 1.8. 

4 Results 

In this section we present the results of our analysis in the MSSM with real and com- 
plex parameters as well as in special scenarios with large mass splittings giving resonant 
enhancement of the signal. 

4.1 MSSM with real parameters 

Unless otherwise noted we use in the following a standard maximal mixing scenario, 
m™ ax , with X t = 2Msusy as defined in [65]. All SUSY parameters needed as input to 
calculate the Higgs masses and mixing matrix with FeynHiggs are given in table 1. 
For the electroweak parameters, m t = 178 GeV, m& = 4.7 GeV, m\y = 80.426 GeV, 
and mz = 91.187 GeV have been used. Gf has been calculated from the running chew 
as Gf = na E \y /(V% sin 2 9 w^w)- We have checked that the impact of using a lower 
m t = 172 GeV on the signal rates and distributions is small compared to the uncertainties 
from higher order effects and can safely be neglected. This also applies to the branching 
ratio Br(H + — > t + u t ), which in the case of m^± = 400 GeV is reduced from 22.7 to 22.0 
%. 

The mass and tan/? dependence of the cross-section after all cuts of table 2 are shown 
in figure 7 as solid curves whereas dashed curves denote the cross-section for the harder 
cuts p± Tjct > 100 GeV and $± > 100 GeV. The left plot is for tan (3 = 50 whereas in 
the right plot we have used m#± = 175 GeV (400 GeV) for the solid (dashed) line. The 
horizontal lines indicate the cross-section needed for -M= = 5. Using this criterion to 
define a detectable signal we see that this would indeed be the case for tan/3 > 30 if 
m H ± = 175 GeV and for 150 GeV < m H ± < 300 GeV if tan (3 = 50 with the softer cuts 




m H ± (GeV) tan (3 

Figure 7: mass and tan (3 dependence of the integrated cross-section in the m™ ax 
scenario. Solid curves are with all cuts of table 2, and dashed curves are with the harder 
cuts p± T - ct > 100 GeV and tf± > 100 GeV. The horizontal lines correspond to -4= = 5. 



p± T - et > 50 GeV and $± > 50 GeV, whereas with the harder cuts tan (3 has to be larger 
than at least 50 in order to have a detectable signal. 

Figure 8 shows the resulting m± distribution for a charged Higgs mass of 175 GeV as 
well as 400 GeV in the case tan/3 = 50 compared to the background after all cuts in 
table 3 have been applied. In the high mass case the harder cuts p± rjot > 100 GeV and 
$± > 100 GeV are used but all the other cuts are the same. For an integrated luminosity 
of 300 fb" 1 and a r detection efficiency of 30% we get S/ \/B = 17 for m#± = 175 GeV and 
S/\fB = 3.2 for m,H± = 400 GeV. (The latter result is slightly better than what would 
have been obtained with the softer cuts p±r jot > 50 GeV and^j_ > 50 GeV.) Here, and in 
all following calculations of S/ \fB, we have used the cut m_|_ > 100 GeV to get S and B. 
In the high mass case S/ \/B could in principle be improved by imposing a harder cut on 
m±, as can be seen from the figure, but not much. Also the possible use of upper cuts 
m± < 200 GeV (m± < 500 GeV) for m H ± = 175 GeV (m H ± = 400 GeV) only marginally 
improves S/y/B from 17 (3.2) to 19 (3.3). In the same figure we also see that the harder 
cuts create a fake peak in the background which could make it more difficult to tell if there 
is a signal. This fake peak appears since the r and v T are mostly produced back to back 
in the W + 2 jets process. Using the harder cuts p± Tict > 100 GeV and fi± > 100 GeV, to 
reduce QCD background and detector misidentifications, the significance for m#± = 175 
GeV and tan (3 = 50 is reduced to S/\/B = 3.1. However, in this case using a upper cut 
m± < 200 GeV is beneficial leading to a significance of S/ \/B = 6.4. This means that 
if the softer cuts p± Tict > 50 GeV and p 1 ^ > 50 GeV are sufficient to suppress the QCD 
background these are clearly to be preferred. 



4.2 MSSM with complex parameters 

In the general MSSM many parameters can be complex. However, our signal process, 
analyzed as described in section 2, is only affected by CP violation in the neutral Higgs 
sector because of the neutral Higgs bosons exchanged in the s-channel. The leading con- 




Figure 8: Comparison of the m_i_ distribution between the signal in the m™ ax scenario 
with tan/3 = 50 and m H ± = 175 GeV (left) as well as m H ± = 400 GeV (right) together 
with the respective backgrounds with all cuts of table 2 (for the high mass the harder 
cuts p±r ict > 100 GeV and/i > 100 GeV are used). The dashed curve denotes the signal 
and the dotted one the background whereas the solid curve is the sum of the two. 



tributions to the CP violation in the neutral Higgs sector arise from loops of the scalar 
top and (to a lesser extend) of the scalar bottom sector where the possibly complex 
Higgs/higgsino mass parameter // and the trilinear scalar couplings A t and are dom- 
inant. Furthermore, in constrained MSSM scenarios implying universality conditions at 
the GUT scale and rotating away all unphysical phases with help of U(l) symmetries of 
the theory, only two phases remain, the phase of ji and a common phase for the trilinear 
couplings [38,39]. Hence, we concentrate in the following on the phases M and (pA t °f A 4 
and A t , respectively, which have the largest effect on the neutral Higgs sector and thus 
possibly affect our signal, assuming <pA b — <Pa t — <f>A t - We have varied M and <pA t inde- 
pendently in the range — n < <fi < ir in order to investigate the phase dependence of our 
signal. However, in the maximal mixing scenario as well as in scenarios with less and no 
mixing in the third generation squark sector, see table 1, we find only small (~ 5%) M 
and (ft A t dependencies of the total cross-section as can be seen in figure 9. 

The small effects of the phases in these MSSM scenarios are basically due to the neutral 
Higgs bosons having quite similar masses {m 2 H _ — m 2 H ^ <C {m H ± + my/) 2 —m 2 Hi ) and small 
widths with the latter also being true for the charged Higgs boson. Bearing in mind that 
s > {m H ± + mw) 2 , the propagators, = l/(s — rn 2 H . + imH^Hi), in Eq. (8) and (9) 
are all of similar size and approximately real. Therefore, as a first approximation, they 
can be put in front of the sum over the different Higgs bosons. The remaining sums of 
coupling factors then reduce to simple numbers, the first to 2 sin 2 (3 and the second to 
—2 sin /3, irrespective of the entries in the Higgs mixing matrix, Oji. Thus, unless we have 
large differences in the masses of the neutral Higgs bosons, the effect from the phases via 
the couplings is small. The phase dependence also enters indirectly via the Higgs masses, 
which typically have a large dependence on the phases. It turns out that in the scenarios 
we are considering the total cross-section has an almost linear dependence on the mass 
of the heaviest neutral Higgs boson and the phase dependence from the masses is much 
larger than that from the couplings. 
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Figure 9: Total cross-section as a function of M and in the maximal mixing scenario 
(left), the less mixing scenario (middle) and the no mixing scenario (right). The light 
shaded areas are in agreement with the constraints from a M and Spo and the dark shaded 
ones are also in agreement with the constraints from EDMs. 



Figure 9 also shows the allowed regions from the following precision measurements, the 
anomalous magnetic moment of the muon, a M , the p-parameter, 5po and the electric dipole 
moments (EDMs) of e, n and Hg. The experimental value of a M from [66] is 

a)f p - = 11 659 208.0 ± 5.8 x 10~ 10 (18) 

and the theoretical value determined using e + e~ data from [67] is 

al M = 11 659 184.1 ±8.0 x 10~ 10 . (19) 

Based on these values we use the following 2a-range for the difference 

-3.7 x 10~ 10 < Aa M < 51.5 x 10~ 10 . (20) 

The experimentally allowed range of 5p is [2] 

-0.0010 < 5p < 0.0025, (21) 

whereas we use the following upper limits on the EDMs [68-70] 

| del < 1-6 • 10~ 27 , (22) 

\d n \ < 6.3 • 10~ 26 , (23) 

\d Hg \ < 2.1 • 10~ 28 . (24) 

In order to calculate a M , Spo, \d e \, \d n \ and I^Hgl in our scenarios and to analyze their 
dependence on the phases of p and A t we have used FeynHiggs 2.2.10, where, however, 
it is important to keep in mind that the EDM routines are not yet fully tested. As can be 
seen from figure 9 the restrictions from the precision measurements are very severe. The 
bound on effectively removes the region with \<f>^\ > n/2 and in the remaining region 
only the light shaded areas are consistent with the Spo bound. Finally the darker shaded 
those consistent also with the EDMs. Hence only small variations from phase 
zero or n are allowed. 
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Table 4: MSSM parameters for the resonant scenario as well as the range for the scan of 
parameters together with the stepsize. 



The constraints from the EDMs are very strong in constrained MSSM scenarios, allowing 
only small values of the phases, especially of M . However, in unconstrained SUSY they 
are rather model dependent, permitting in general larger phases. For example, due to 
cancellations between different SUSY contributions to the EDMs or in SUSY models with 
heavy sfermions in the first two generations larger phases may be allowed [71-79] . Recently 
it has been pointed out that for large trilinear scalar couplings A, phases 0^ ~ 0(1) can 
be compatible with the bounds on d e , d n and d^ g [80]. Furthermore the restrictions on 
the phases may also disappear if lepton flavour violating terms in the MSSM Lagrangian 
are included [81]. In conclusion this means that large phases cannot be ruled out and 
therefore we analyze the full range of the phases to determine possible effects. 

To see the dependence of the allowed regions on the chosen SUSY scenario we show in 
figure 9 also the results for the less mixing and no mixing scenarios. In the less mixing 
scenario the constraints from 5po are less severe but the EDMs still give hard constraints 
whereas for the no mixing scenario we get a complete band around M = 0. Finally we 
have also checked that similar results are obtained for tan (3 = 10. 

In the MSSM scenarios studied in this section the CP-odd rate asymmetry, equation (13), 
is always quite small, \A C p\ < 0.3%. This due to the fact that in order to get a large 
asymmetry the effects from the absorptive phases in the Higgs propagators, S Hi as well 
as the phases in the couplings both have to be large. On the contrary, in the scenarios 
considered here the phases in are always quite small and at the same time the mixing 
between the neutral Higgs bosons is also typically small resulting in small asymmetries. 
In order to get a non- negligible asymmetry one needs large phases from at least one of the 
propagators and at the same time large mixing between the neutral Higgs bosons, which 
seems difficult to achieve in the MSSM. In more general 2HDMs, large asymmetries should 
be possible in scenarios with two resonant neutral Higgs bosons in the s-channel and large 
mixing, although this remains to be verified. 

4.3 Resonant scenarios 

In SUSY parameter regions with \A t \, or \A b \ > 4Msusy the dominant terms of the 
1-loop corrections to the quartic couplings in the Higgs sector [39] can induce a large 
mass splitting between the charged and neutral Higgs bosons [82,83]. For example, in the 
scenario given in table 4, the CP-odd Higgs is 80 GeV heavier than the charged Higgs 
allowing resonant production in the s-channel. 

In order to examine the perturbative stability of this scenario with very large 1-loop 



corrections to the mass of the CP-odd Higgs boson we list in table 5 the masses of the 
Higgs bosons calculated with FeynHiggs 4 at tree-level, at 1-loop order, and with all 
available corrections. For comparison, the Higgs masses in the maximal mixing scenario 
are also given. From the table it is clear that the higher order corrections have much 
smaller impact than the leading ones suggesting that the perturbative expansion is under 
control. We have also verified that similarly large mass splittings are obtained with 
CPsuperH using all available corrections although not in precisely the same parameter 
points. The latter feature is due to this kind of resonant scenarios being rather fine-tuned 
and therefore sensitive to differences in the implementation and approximations used in 
the two programs. 
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Table 5: Masses of the Higgs bosons in the maximal mixing and resonant scenarios, cal- 
culated with FeynHiggs at tree-level, at 1-loop order and using all available corrections. 

To get an indication of how fine-tuned this kind of resonant scenario is, we have performed 
a scan over the relevant parameters as shown in table 4. Defining a resonant scenario via 
the relation > m H ± + mw, there is in fact a large range in both tan (3 and m#± 
were such scenarios are found as can be seen in figure 10. In each tan/3 and m^± point 
about 9300 different scenarios are tested. On average about one half of these correspond 
to physical scenarios and in turn about 0.2% of the latter give resonant conditions, which 
illustrates the level of fine-tuning in these scenarios. The dependence on tan (3 and m#± 
is illustrated in figure 10. At the same time, these scenarios typically have relatively low 
squark masses so for large m#± (> 200 GeV) the decay to squarks becomes dominant and 
thus the specific analysis we present here is not suitable. At the same time, in scenarios 
with light squarks there may also be a resonant enhancement in the gluon initiated channel 
from squark loops as already discussed above. This occurs if the sum of the squark masses 
are close to threshold. In addition is is also possible in this case to have enhancement 
from the s-channel resonance. Thus in the rare situation that both these effects occur 
simultaneously, without opening the decay channel of the charged Higgs into squarks, the 
signal could in fact be enhanced. Finally, we note that the study of resonant scenarios 
also illustrates what could happen in a general 2HDM, where the masses are more or less 
independent parameters. 

In case of resonant production, the H ± and W bosons are produced with typically small 
transverse momenta. Thus it is favourable to loosen the cuts on the light jets from the 
W . Applying the basic and additional cuts from table 2, except the cuts p±hj > 50 GeV 
and p± S j > 25 GeV on the light jets, we get an integrated cross-section of 52 fb for a 
charged Higgs boson mass of 175 GeV in the resonant scenario given in table 4. Figure 11 
shows the resulting m_i_-distribution compared to the background. With an integrated 
luminosity of 300 fb -1 and a r detection efficiency of 30% we get a significance S/y/B = 

4 We have added the partial decay width T(Hi — > i^W^) to the calculation of T Hi in FeynHiggs. 
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Figure 10: Contour plot of the relative number of resonant scenarios to physical scenarios 
in each point. The lines are 0.01%, 0.05%, 0.2% and 1%. The region around tan/3 8 
contains more physical scenarios than the upper region leading to a total of 0.2% resonant 
scenarios. 

56. For comparison, if we apply the harder cuts p± Tjct > 100 GeV and $±_ > 100 GeV 
in this resonant case the significance is reduced drastically to S/VB = 0.2 due to the 
typically small transverse momentum of the if ± -boson. Thus, in the case of harder cuts 
the resonantly enhanced cross-section is only of use if m#± is large enough such that 
m H ±/2 is well above 100 GeV. 

Finally we have also tried to investigate the phase dependence of the cross-section in this 
resonant scenario. Figure 12 shows the dependence of the cross-section on </> M , where the 
Higgs masses, couplings and widths have been calculated with FeynHiggs at one-loop 
accuracy 5 . Note, that the phase 4>A t is irrelevant in this scenario because \A t \ = 0. For 
comparison, the figure also shows the result with all available corrections in the case 
M = 0. The very large phase dependence is due to the fact that the production goes 
from non-resonant to resonant when varying the phase. More specifically, as can be seen 
in table 5 we get m A = m H . } = 246 GeV for M = in the 1-loop case, which is below the 
resonant threshold, whereas m# 3 = 342 GeV for the largest values of M where we got a 
stable result, which is clearly in the resonant regime. Even in these resonant cases, the 
width of the H% is typically small, Y/m < 0.02, and consequently the CP asymmetry is also 
small. However, we have not been able to make a more thorough study of the sensitivity 
to these phases due to numerical instabilities nor have we searched for resonant scenarios 
with large Y/m as well as large mixings since this requires a dedicated study. 

5 Summary and conclusions 

In this paper we have studied the viability of detecting charged Higgs bosons produced 
in association with W bosons at the upcoming LHC experiments. Since our study is 
of exploratory character we have stayed on the parton level, although with appropriate 



5 A calculation with all available corrections is not possible here because the phases lead to numerical 
instabilities. 
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Figure 11: Comparison of m± distribution between signal (for the resonant scenario with 
m H ± = 175 GeV and tan/3 = 11 given in table 4) and background with all cuts of table 2 
except p± h j > 50 GeV and p± sj > 25 GeV. The dashed (dotted) curve denotes signal 
(background) and the solid curve is signal and background combined. 
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Figure 12: Total cross-section, with the Higgs masses etc. calculated to one-loop, as a 
function of M for the resonant scenario. For comparison the result with all available 
corrections for M = is also shown (labeled 2- loop). The light shaded areas are in 
agreement with the constraints from a M and Spo and the dark shaded ones are also in 
agreement with the constraints from EDMs. In the black areas numerical instabilities 
occurred in the calculation of the Higgs masses and mixing matrix. 



smearing of momenta, and leave the inclusion of parton showering, hadronisation, full 
simulation of the detector, etc. for future studies. On the same vain we have not tried to 
include any higher order corrections in the production cross-section. The only exception 
is that we use running quark masses in the Yukawa couplings, which gives a cross-section 
that is almost a factor three smaller compared to when using the pole-mass, in agreement 
with NLO QCD [31]. Other higher order corrections on the production cross-section, such 
as the choice of factorisation scale, are expected to be O(10 — 20%). 

As we have shown in this paper, using the leptonic decay of the charged Higgs boson 
and hadronic W decay, giving rise to the signature Tj et $± + 2 jets, it is possible to design 
appropriate cuts against the irreducible Standard Model background from W+ 2 jets 
production. This is in contrast to earlier studies using hadronic decays of the charged 
Higgs bosons, which was found to suffer from too large irreducible background due to tt 
production [26]. At the same time we find that the significance of the resulting signal 
does depend on the assumptions made on the cuts needed against reducible backgrounds 
(mainly from QCD and detector misidentifications). 

In the standard maximal mixing scenario of the MSSM, using the softer cuts (p± r t > 
50 GeV and $ ± > 50 GeV) we find a viable signal in the case of large tan (3 (> 30) 
and intermediate charged Higgs masses (150 GeV < ttlh± < 300 GeV), whereas with 
the harder cuts (p± Tjet > 100 GeV and p'± > 100 GeV) we only find a viable signal for 
tan/3 > 50 if in addition an appropriate upper cut on m± = ^2p± r tf±[l — cos(A</>)] is 
applied. Thus, in the best case the associated charged Higgs and W boson production 
could serve as a complement to production in association with top quarks in the difficult 
transition region m#± ~ m t although only for large tan/5, but before one can draw any 
firm conclusions the effects of reducible backgrounds have to be studied in more detail. 

In MSSM scenarios with less or no mixing in the third generation squark sector we get 
similar results as in the standard maximal mixing scenario. In fact, the differences com- 
pared to the maximal mixing scenario are in both cases smaller than the O(10 — 20%) 
effects we expect from higher order electroweak and QCD corrections as well as other 
variations of the SUSY scenario. Similarly we have also found that the cross-section de- 
pends only weakly on the CP-violating phases of the SUSY parameters, even those of A t 
and /i, which result in a mixing between the CP-even and CP-odd Higgs bosons. This is 
a general feature that will be true as long as the differences in the neutral Higgs bosons 
masses are small compared to the charged Higgs plus W-boson mass (more specifically 
m 2 H _ —f^Uj ^ ( m H± + m w) 2 ~ m ' 2 H i ) even if the Higgs mixing matrix is highly non-diagonal. 
The phases of A t and ji also lead to differences in the masses of the neutral Higgs bosons. 
In fact, we have seen that the main effect on the cross-section comes from these kinematic 
effects of the Higgs masses and not from the changes in the couplings due to the phases. 

We have also studied a class of special resonant scenarios where > m#± + m\y for 
one of the neutral Higgs bosons (the CP-odd Higgs A in the MSSM with real parameters) 
leading to resonant production in the s-channel. This results in a very large enhancement 
of the total cross-section (up to a factor 100) in the region of intermediate tan/3 (~ 10). 
However, due to the different kinematics of resonant production, the significance of the 
signal depends very strongly on the p± Tjct and p"± cuts. With the softer cuts we get a 
significance of order 50 in the case of m H ± = 175 GeV, whereas with the harder ones 
it drops to 0.2. Another problem with these type of MSSM scenarios is that for larger 
charged Higgs masses (m H ± > 200 GeV) the decays to squarks opens up which would 



require a different type of analysis. 

One may also worry about the perturbative stability of these resonant scenarios with very 
large one-loop corrections to the mass of the CP-odd Higgs boson. However, comparing 
the one-loop and available two-loop corrections, as was done in Table 5, this does not 
seem to be the case. When making a sparse scan of parameters we find similar scenarios 
in a large range of tan/5 and m H ± and we also get similar results when using CPsuperH 
instead of FeynHiggs although not precisely at the same points in parameter space. 
We have also found that in these resonant scenarios the cross-section can have a large 
dependence on the CP violating phases of the SUSY parameters but we have not found any 
appreciable CP asymmetries. However, we have not been able to make a more thorough 
study of the sensitivity to these phases which would require a separate study. 

Finally, it should be emphasized that our study is specifically for different MSSM scenarios 
and that the conclusions may change in other models. For example, in a general 2HDM, 
the resonant enhancement is more natural since the Higgs masses are independent of 
each other. In addition there will be no charged Higgs decays to squarks as we found 
in the MSSM thus making it possible to have a clear signal for large charged Higgs 
masses even with the harder cuts against reducible backgrounds. Similarly nonminimal 
supersymmetric models with larger Higgs sectors may also offer more natural possibilities 
for resonant enhancement. 
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